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Lattice results in supersymmetry are summarized. Past, present and future perspectives are discussed. 



1. INTRODUCTION 

Supersymmetry or fermion-boson symmetry is 
one of the most fascinating topics in field theory. 
Even if has not yet been observed in Nature, thou- 
sands of papers have been written on this subject. 
Its validity in particle physics follows from the 
common belief in unification through the feasi- 
bility of incorporating quantum gravity. From a 
theoretical point of view, non-perturbative stud- 
ies of supersymmetric (SUSY) gauge theories turn 
out to have remarkably rich properties which are 
of great physical interest, as have been shown in 
|l|. For this reason, much effort has been ded- 
icated to formulating a lattice version of SUSY 
theories (for a recent review of SUSY Yang-Mills 
(SYM) theories, using Wilson fermions, see Q). 
The lattice formulation has been succesful to ex- 
tract non-perturbative dynamics in field theory, 
specially in QCD, and may be able to provide 
additional information and confirm the existing 
analytical calculations. Whether superymmetry 
is or not an exact symmetry is a question that 
must be settle by going beyond perturbation the- 
ory. 

1.1. SYM theory in the continuum 

In the continuum, the action for the N = 1 
SYM theory with an SU (N c ) gauge group is 



£ = -\f; v {x)f« v {x) + ±\ a ( x ) 1 ^\( x y , (i) 



where A a is a 4-component Majorana spinor and 
satisfies the Majorana condition A a = X aT C. The 
gluon fields are represented by A p — ~igA p T a 
and F pu = -igF^T a . V p X a = d p X a + gf abc A^\ c 



is the covariant derivative in the adjoint repre- 
sentation. The continuum SUSY transformations 
read 



SAJx) = -2g\(x)j„e, 



(2) 



6\{x) 



-a pT F pT (x)e , 8\{x) = -ea pT F pT (x) , 



where a pT = |[7p,7r]j A = \ a T a and e is 
a global Grassmann parameter with Majorana 
properties. These transformations relate fermions 
and bosons, leave the action invariant and com- 
mute with the gauge transformations so that the 
resulting Noether current S p (x) is gauge invari- 
ant. For N = 1 SYM theory the supercurrent is 
given by 



S li (x)=-F a (x)a prlli \ a {x). 



(3) 



Classically, the Noether current is conserved, 
d p S p (x) = 0, provided the fields satisfy the equa- 
tions of motion. Furthermore one finds that it 
fullfills a spin 3/2 constraint: r y p S p (x) = 0. 

1.2. Superfields 

Supersymmetry in Eq. (|l|) can be better un- 
derstood in the language of superfields . In 
the Wess-Zumino gauge, the action for the N = 1 
SYM theory is 



S SYM = Re <^ / d 4 x d 2 9 Tr[W a W a 
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(4) 



W(x,6,6) a = -|(D D)e~ 2V D a e 2V is the spino- 
rial field strength superfield which depends on 
the four coordinates x and the anticommuting 
Weyl spinor variables 9 a , 9 a (a, a — 1,2). D a 
and D a are the superspace derivatives while V = 
V(x, 9, 9) is the vector superfield which takes the 
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form j| 

-i(ee)0^ w (x) + ^(90){00)D{x) . (5) 

The particle content can be identified as one 
bosonic vector field A^ = -igA^T a , its SUSY 
partner, a complex Weyl spinor field i/) w = T" 1 ^, 
and an auxiliary scalar field D. The action can 
be expressed in terms of these components as 

ssym = J d 4 z j - l -F^ u F^ a + \ir w ^{v^ w r 
-^(^r^c + ^ 2 |- (6) 

In the on-shell case, the auxiliary scalar field D 
is eliminated by the equation of motion D = 
0. Introducing one massless Majorana fermion, 
the gluino, in the adjoint representation A = 
((jpw)a, (ipw) a )t and going to the Euclidean space, 
we recover Eq. ([!]). 

1.3. SUSY Ward Identities (Wis) 

The existence of a renormalized SUSY Ward- 
type identity 

dpSfc) = 2m RX B.(x) , (7) 

is generally assumed. It is obtained by perform- 
ing a variation of the functional integral with re- 
spect to a local, smooth transformation, e = e(x), 
and then putting the sources to zero. is the 
renormalized supercurrent and \R = Z x x, with 
X = ^F"o ' M i/A° ' . mR is the renormalized gluino 
mass. For niR — we have supersymmetry while 
a non- vanishing value of m r breaks supersymme- 
try softly. 

1.4. Non-perturbative effects of the SYM 
dynamics 

Because of asymptotic freedom and the corre- 
sponding infrared instability, it is reasonable to 
guess that, as in conventional QCD, there are 
many non-perturbative phenomena that occur in 
this theory. For the N ~ 1 SYM theory we may 
expect confinement and spontaneous chiral sym- 
metry breaking. The discrete chiral symmetry is 



expected to be broken by a non-zero gluino con- 
densate while the confinement is realized by color- 
less bound states described by the effective action 
belonging to chiral supermultiplets. 

It is generally assumed that supersymmetry is 
not anomalous (Eq. (Q) holds) and only the mass 
term is responsible for a soft breaking. However, 
in Q the question of whether non-perturbative 
effects may cause a supersymmetry anomaly has 
been raised. Only a study of the continuum limit 
of the lattice SUSY Wis can shed light on this 
question. 

1.5. Chiral symmetry breaking 

Introducing a non-zero gluino mass term in 
Eq. (|l|), Cmass = nig\ a \ a , breaks super- 
symmetry softly (which implies that the non- 
renormalization theorem and cancellation of di- 
vergencies are preserved Q). In the massless case, 
the global chiral symmetry is U(l)\: 

A^e-^A, A — > Ae^ 75 , (8) 

which is anomalous. This implies that the diver- 
gence of the axial current, = A75A, is 

f) J 5 — ~^ cg r V-vpa p« J?a (q\ 

m m ~ 32tt 2 r ^ r P a- \v) 

The anomaly leaves a Z 2 n c subgroup of £7(1) a un- 
broken. Eq. (^) is equivalent to rrig — » m,ge~ 2l(p ' t5 
and Qsym -> ®sym - 2N c tp. 

In the SUSY case, m g = 0, the U(l)\ symme- 
try is unbroken if <p = ^- for (k = 0, 1, • • • , 2N C - 
1)- Z 2 n c is expected to be spontaneously broken 
to Z 2 by a value of ( AA) ^ f) . The consequence 
of this spontaneous chiral symmetry breaking is 
the existence of a first order phase transition at 
rrig = 0. That means the existence of N c degen- 
erate ground states with different orientations of 
the gluino condensate (k = 0, • • • , JV C — 1), 

(AA) = cA 3 e^ , (10) 

where A is the dynamical scale of the theory 
which can be calculated on the lattice, for ex- 
ample, while c is a numerical constant which 
depends on the renormalization scheme used to 
compute A. Eq. ( |l(i| ) shows the dependence on the 
gauge group. For SU{2) two degenerate ground 
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states with opposite sign of the gluino conden- 
sate, (AA) < and (AA) > appear, while for 
SU(3) there are three degenerate vacua at k = k c 
(in g a fourth state with (AA) = is claimed). A 
first numerical study of SYM theory with a gauge 
group SU(3) is in §. 

1.6. Magnitude of the gluino condensate 

The calculation of the gluino condensate for 
the N = 1 SYM theory is a puzzle. Two ap- 
proachs in the literature have been used which 
give different results for c. One is based on 
weak-coupling instanton (WCI) calculations [|l0| 
and gives c = 1. In the second, calculations 
based on strong-coupling instanton (SCI) give 
c = 2/((N c - l)!(3iV c - l)) 1 /^ [JlTJ] . For a review 
on SCI see [|l2|. Several suggestions have been 
put forward to resolve the puzzle, see for exam- 
ple [|l3j. Ref. |l4|] cast serious doubts on the SCI 
calculation by showing that the cluster decompo- 
sition is not valid. It is also shown that the addi- 
tion of the so-called Kovner-Shifman (KS) chiral 
symmetric vacuum state ||[l5) can not straight- 
forwardly resolve the disagreement between the 
SCI and WCI results @. The KS vacuum can 
indeed potentially resolve the mismatch at the 1- 
instanton sector (k = 1) but it fails to do so for 
the topological sectors with k > 1 Jm). Using an 
instanton liquid picture gives qualitatively simi- 
lar results and evidence for the gluino condensate 
0. 

1.7. Light hadron spectrum 

For the N — 1 SYM theory a low energy effec- 
tive action has been proposed by Veneziano and 
Yankielowicz (VY) |L7|. The action contains all 
degrees of freedom, gauge invariant and colorless, 
composite fields: F« V F* V , F* V F* V , A a A a , in anal- 
ogy to QCD, while \ = a ^F^ v \ a is a new type 
of composite operator formed by a gluino A and a 
gauge field F, both in the adjoint representation. 
These fields can be combined to form the chiral 
supermultiplct, containing the expression for the 
anomalies as component fields p3 



S(x, 9) = <p(x) + V2d X (x) + 69F(x) . 



(11) 



4> = @2g('4>w) '(' t ( ) w)aj is proportional to the 
gluino bilinear (here ip w denotes a 2-component 



Weyl spinor), while the other components, which 
we do not report here (see J3|), contain combina- 
tions of gluino-gluino field and gluino-gluon fields. 
The VY action has the form [|| 

C eff = 1(5^)1/3^ + 7 [(Slog4r ~ S)\ F + h.C.] . 

a /i J 

Expanding the effective VY action around its 
minimum, it is found that the low energy spec- 
trum forms a supermultiplet, consisting of a 
scalar meson A a A a (the a — /o), a pseudoscalar 
meson A a 75A a (the a — r]/), (the a denoting the 
adjoint representation), and a spin- 1/2 gluino- 
glueball particle, the x- Glueballs are absent in 
this formulation. 

In the SUSY point these masses are degener- 
ate. The introduction of a nig ^ breaks super- 
symmetry softly and leads to a splitting of the 
multiplet. How the spectrum is influenced by the 



soft SUSY breaking has been studied in [19 
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3Nr. 



(12) 



Unfortunately, the range of applicability of the 
linear mass formulae is not known because of the 
unknown magnitude of the constants and of the 
higher order terms in Eq. (^2[). 

Introducing an extra term in the effective VY 
action, Farrar et al. (FGS) solved the question 
of including glueballs in the low energy spectrum 
PQ| . For unbroken supersymmetry we expect to 
see two chiral multiplets (not one), at the bot- 
tom of the SYM spectrum. The lighter one con- 
tains a + glueball « F^F^, a 0~glueball w 
e iivpo 'F^Fpcr and a gluino-glueball ground state, 
while the heavier supermultiplet is the VY one. 
Moreover, in the FGS picture, a non-zero mixing 
between the a — fo and + glueball is possible. 
Of course, other generalizations beside the FGS 
picture are conceivable. 
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2. LATTICE FORMULATION OF N = 1 
SYM 

The question of whether it is possible to for- 
mulate successfully SUSY theories on the lattice 
has been addressed in the past by several au- 
thors |2l],^,^3| (the reader is refered to Q for a 
sucessful contruction of a lattice 2c? Wess-Zumino 
model, and some discussion concerning Ref. p2[). 
It can be seen that a lattice regularized version 
of a gauge theory is not SUSY since the Poincare 
invariance (a sector of the superalgebra) is lost, 
thus {QonQp} = ^CafiPfj,- Poincare invariance is 
achieved automatically without fine tuning in the 
continuum limit because operators that violate 
Poincare invariance are all irrelevant. Moreover, 
if the SUSY theory contains scalar fields one can 
have scalar mass terms that break supersymme- 
try. Since these operators are relevant, fine tuning 
is necessary in order to cancel their contributions. 

Another problem is the question of how to bal- 
ance bosonic and fermionic modes, the numbers 
of which are constrained by the supersymmetry: 
the naive lattice fermion formulation results in 
the doubling problem p4| , and produces a wrong 
number of fermions. The problem can be treated 
as in QCD by using different fermion formula- 
tions. Let us briefly summarize those which have 
applications in SUSY theories. 

2.1. Wilson fermions 

In the Wilson formulation for the N = 1 SYM 
theory it is proposed to give up manifest super- 
symmetry on the lattice and restore it in the con- 
tinuum limit Supersymmetry is broken by 
the lattice, by the Wilson term and is softly bro- 
ken by the presence of the gluino mass. Super- 
symmetry is recovered in the continuum limit by 
tuning the bare parameters g and the gluino mass 
nig to the SUSY point. The chiral and SUSY 
limit can be recovered simultaneously at m§ = 0. 

In the Wilson formulation, the Curci and 
Veneziano (CV) effective action, suitable for 
Monte Carlo simulations, reads 

Scv = f 1 - TT ReTrC v) - ^logdetQ[[/] . 
For the gauge group SU(N C ), the bare coupling 



is given by (3 = 2N c /g 2 . The fermion matrix for 
the gluino Q is defined by 



(yb 



4 

,xa[U] = 5 xy S ab -ky][ 



(13) 



$y,x+p,(l + 7m)^<w + <W,*( 1 ~ Tm)^,^] • 

k is the hopping parameter defined as k = 
1/(2(4 + mofl)), where tuq is the bare mass, and 
the matrix for the gauge field link in the adjoint 
representation is 



[U] = lTr(U^T a U Xfl T b ) . (14) 



The fermion matrix for the gluino in Eq. (|13|) is 
not hermitian but it satisfies — 75Q75. That 
allows for the definition of the hermitian fermion 
matrix Q = 75 Q. The path integral over the Ma- 
jorana fermions gives the Pfaffian, 

[d\] e -? XQX = f [d\]e-? xTcQX = Pf(M) , 



where M = CQ is an antisymmetric matrix. 



It is easy to see that Pf(M) = ±y/det Q. In 
the effective CV action the absolute value of the 
Pfaffian is taken into account (this may cause the 
sign problem). The omitted sign can be included 
by reweighting the expectation values according 
to the formula, 



(0) = 



(OsignPf(M)) 



cv 



(signP/(M)) 



(15) 



cv 



The spectral flow is a method which checks the 
value of the sign of the Pfaffian. The experience of 
the DESY-Munster collaboration shows that bel- 
low the critical line fc c (/3), corresponding to zero 
gluino mass (rn§ = 0), negative Pfaffians practi- 
cally never appear |32| , ^3| , p3| . 

The factor 1/2 in front of logdetQ[J7] shows 
that Majorana fermions imply a flavor number 
Nf — 1/2. This can be achieved by the hy- 
brid molecular dynamics (HMD) algorithm J2j| 
which is applicable to any number of flavors. 
The HMD has been checked for the CV action 
in N = 1 SU{2) SYM at small lattices 4 3 x 8 
p6| |. Another method for simulation with non- 
even numbers of flavors is based on the multi- 
bosonic algorithm proposed by Liischer p7|. A 
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two-step variant using a noisy correction step p8| , 
has been developed by Montvay in |2^,0 called 
the two-step multibosonic (TSMB) algorithm. In 
the two-step variant, to represent the fermion 
determinant one uses a first polynomial V^}{x) 
for a crude approximation realizing a fine correc- 

(21 

tion by another polynomial Vn 2 ( x ) that satisfies, 
lim V$(x)V®(x) = x~ N f/ 2 , for x € [e,A]. 
The fermion determinant is approximated as [H 

1 



10 _ lattice: 6 3 » 12 



det(QtQ) 



N f 



detP^ ) (Qtg)detP^ ) (QtQ) 



Unquenched results using TSMB have been re- 
ported by the DESY-Mimster-Roma collabora- 
tion (the first large scale numerical simulation of 
N = 1 SYM theory). For SU(2) see §|||||, 
while for SU(3) see |J. Interesting quenched re- 
sults (which was pioneering work) are in p^ . p5[ 

D- 

In order to check the expected pattern of spon- 
taneous chiral symmetry breaking, let us first 
write the expression for the renormalized gluino 
condensate, obtained by additive and multiplica- 
tive renormalizations: (AA)jj/^ = Z(a(i)[(\\) — 
bo(afi)]. In a numerical simulation, a first or- 
der phase transition (or cross-over) should show 
up (on small or modere lattices) as a jump in 
the expectation value of the gluino condensate at 
k = k c . By tuning the hopping parameter k to 
k c , for a fixed j3 value, one expects to see a two 
peak structure in the distribution of the gluino 
condensate. By increasing the volume the tun- 
neling between the two ground states becomes 
less and less probable and at some point prac- 
tically impossible. Outside the phase transition 
region, the observed distribution can be fitted by 
a single Gaussian but in the transition region a 
good fit can only be obtained with two Gaussians. 
The DESY-Miinster collaboration performed the 
first lattice investigation of the gluino conden- 
sate. Results for SU(2) are reported in Fig. [l] 
and show that a first order phase transition oc- 
curs at k c = (0.1955 ± 0.0005) for /3 = 2.3. The 
lattice size used is 6 3 x 12 pi| . For SU(3), prelim- 
inary results are in [^) . Here the pattern is more 
complicated, Z§ — > Z2. The probability distribu- 
tion for the gluino condensate was measured for 
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Figure 1. Distribution of the gluino condensate 
for k = 0.195 and k = 0.196. k c = (0.1955 ± 
0.0005). 



k = 0.1950 and (3 = 5.6 on a rather small lat- 
tice 4 3 x 8, with encouraging results. Of course, 
a more accurate study by going to larger lattices 
will clarify the nature of the transition. 

SUSY restoration can be also verified by a di- 
rect inspection of the low energy mass spectrum 
]32j , |42] ,Pl : this is expected to reproduce the SUSY 
multiplets predicted in . An accurate anal- 

ysis of the spectrum is, however, a non-trivial task 
from the computational point of view and an in- 
dependent method for checking SUSY restoration 
would be welcome. Moreover, the mixing between 
the states has been measured. Numerical simu- 
lations show that this mixing is practically zero 
|3^,|3£|. At this conference, new results for the 
low lying spectrum of N = 1 SYM, for 12 3 x 24 
and k = 0.1955 and for 16 3 x 32 and k = 0.194, 
are reported by Peetz |yj . 

A study with a large number of colors N c and 
strong coupling, by considering the hopping pa- 
rameter expansion as a sum over lattice paths 
(random walks), is |p^| . The surprising result is 
that the quenched approximation is exact at this 
order (and it is also consistent with the quenched 
results of Donini et al. in |35|j ). Formulae for 
the propagators and masses of 2- and 3-gluino 
states are also presented. The 2-gluino masses 
do coincide with the results for the meson spec- 
trum in ordinary lattice QCD at strong coupling 
]37| ]. A preliminary study of a 3-gluino states is 
reported also in . This particle does not ap- 

pear in p7|,pQ| even if it contains the same quan- 
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turn numbers. A numerical analysis of this issue 
would clarify whether it contributes or not in the 
mass spectrum. At this conference, a study of the 
strong coupling expansion for N — 1 SYM theory 
using the Hamiltonian formalism is also presented 
pq| . It is shown that the theory is effectively de- 
scribed by an 5*0(4) antiferromagnet. 

2.2. Domain wall fermions 

A very nice innovation is the domain wall 
fermions (DWF) approach. A new lattice fermion 
regulator, which improves the lattice formula- 
tion for fermions because the zero gluino mass 
is achieved without fine tuning. Application of 
DWF in SUSY theories has been explored in 
§HQ and also suggested in [fH, with a differ- 
ent approach as [Q. First Monte Carlo simu- 
lations for JV = 1 SUJ2) SYM with DWF, using 
the lines of Re fs. B fl are in DWF were 

introduced in |48|| and were further developed in 
@@. For a recent review on DWF for SUSY 
gauge theories see |3l) . 

There are two unwelcome difficulties in using 
Wilson fermions as mentioned in the previous sec- 
tion. The first one is the need for fine tuning. The 
second one is related with the sign of the Pfaffian. 
DWF are defined by extending the space-time to 
five dimensions. Also a non-zero five dimensional 
mass or domain wall height mo, which controls 
the number of flavors, is present. L s is the size of 
the fifth dimension and free boundary conditions 
for the fermions are implemented. As a result 
the two chiral components of the Dirac fermion 
are separated with one chirality bound exponen- 
tially on one wall and the other on the opposite 
wall. For any value of a the two chiralities mix 
only by an amount that decreases exponentially 
as L s — > co. For L s = co, chiral symmetry is ex- 
pected to be exact even at finite lattice spacing. 
So, there is no need for fine tuning f5lj . DWF 
offer for the first time the oportunity to separate 
the continuum limit, a — > 0, from the chiral limit, 
L s — > co. 

DWF introduce two extra parameters: L s and 
mo . These two parameters together with the four 
dimensional mass m / control the effective fermion 
mass m e f / . In the free theory one finds ]52| ] 

m eff = ™o(2 - m )[m/ + (1 - m ) Ls ] , (16) 



\\ 8*. 0=2.3, 0.00432(22) 



^ ---- 



4 , 0=2.1, 0.00268(19 



20.0 30.0 
L, 



Figure 2. Dynamical gluino condensate at m/ 
vs L s on two different lattices. 



with < mo < 2. The value of mo = 1 is op- 
timal because finite L s effects do not contribute 
to m e f f . In the interactive theory, one would not 
expect such optimal value, due to the fact that 
mo will fluctuate. Then the goal would be to 
have L s large enough to have the second term of 
Eq. ( [llf ) small, in order to simulate at reasonably 
small masses and extrapolate to the chiral limit, 
m/ — > and L s — > co Q . 

The effective action, which is not reported here 
(more details in jl6|) contains L s heavy species. 
When L s is increased they may introduce bulk 
effects which must be substracted. This can be 
done by introducing Pauli-Villars type fields Q . 
The fermionic path integral gives the Pfaffian 



[d\]t 



-S F 



Pf(M F ) = y/det(D F ) 



(17) 



which is positive for mf > 0. This is a second 
advantage of the DWF approach. 

Numerical simulations with DWF for N = 1 
SU(2) SYM are reported in («|[47|, using the 
(HMDR) of H. They concern only the study of 
the gluino condensate, for rather small volumes, 
4 4 and 8 4 . The results for m/ = are shown 
in Fig. ||. As can be seen, (AA)(L S — > co) has a 
non-zero VEV which remains different from zero, 
even for rather small volumes. These results (and 
also those in [Q) support the non-zero value for 
the gluino condensate. 

Also the DWF has difficulties in its implemen- 
tation. Beside the two extra parameters, it may 
happen that the two chiralities do not decouple, 
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even for L a — > oo. In this case the chiral sym- 
metry can not be restored. This may need large 
values of i s and for this reason much expensive 
cost in the simulations. From the computational 
point of view, SYM is harder to simulate than 
QCD, using DWF, while with Wilson fermions, 
SYM is easier to simulate than QCD. 

2.3. SUSY Wis on the lattice 

Another independent way to study the SUSY 
limit is by means of the SUSY Wis. On the lat- 
tice they contains explicit SUSY breaking terms 
and the SUSY limit is defined to be the point in 
parameter space where these breaking terms van- 
ish and the SUSY Wis take their continuum form. 
These issues have been investigated numerically 
by the DESY-Munstcr-Roma collaboration, using 
Wilson fermions p^ , p3[ . Previoulsy in |3(|. 

In the continuum, the SUSY Wis are given by 
Eq. (0) where S*(x) = Z s S^x) + Z T T^{x) and 
Zs and Zt are multiplicative renormalization fac- 
tors. Tfj,(x) = F^ v (x)'y / j,X a (x) is a mixing current 
with dimension 7/2 as S f _ l (x) fj3f . 

In the Wilson formulation, supersymmetry is 
not realized on the lattice. One might still define 
some lattice SUSY transformations (which reduce 
to the continuum ones (^) in the limit a — > 0). 
One choice is J|3|,|5| 

SU^x) = -agJ7 M (x)e(a;)7 Al A(a;) 

—age(x + a//)7 M A(x -I- afijU^x) , 



SX(x) = --<7 pT Qpr{x)e{x) . 



(18) 



In the case of a gauge invariant operator insertion 
0{y), we find for the bare SUSY Wis 

(0(y)VMx)) - 2m (O(y) X (x)) + (S^|e=o) 

= (O(y)Xs(x)) . (19) 

Comparing with Eq. (Q), beside the presence of 
a non-zero bare mass term in the action which 
breaks supersymmetry soflty, the rest of the 
SUSY breaking results in the presence of the X$ 
term, that appear since the action is not fully 
SUSY. In order to renormalize the SUSY Wis a 
possible operator mixing has to be taken into ac- 
count. In the case of gauge invariant operator 



insertion, X$ mixes with the following operators 
of equal or lower dimension ^4j, V^S^, V M T M and 
X- The SUSY Wis can be written as @ 

= m R Z s 1 {0(y) X (x)) . (20) 

Here the gauge invariant operator 0{y) at point 
y is assumed to be sufficiently far away from x in 
such a way that contact terms are avoided. This 



is the on-shell regime, mji 
V M is the lattice derivative. 



mo 
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In numerical simula- 
tions, Eq. (20) can be computed at fixed j3 and 
k. Thus, by choosing two elements of the 4x4 
matrices, a system of two equations can be solved 
for ZtZq 1 and mnZg 1 . One must clearly ensure 
that these two equations are non-trivial and in- 
dependent. To do this, a proper definition of the 
supercurrent operator and the operator insertion 
O is necessary. In Ref. J33|, two different defini- 
tions for the supercurrent operator are considered 
in Eq. (|20|): the local definition (||), and more 
involved, the point-split definition, reported also 
in (55), which differs from the local one in order 
0(a). 

The dependence of ZtZ s 1 and muZ s 1 on k, 
are shown in Fig. || and Fig. [|, respectively, with 
two different operator insertion, x(x) and Tq(x), 
for a lattice size 12 3 x 24. Fitting results are con- 
sidered only for the insertion x( x )- hi Fig. 
ZtZ^ 1 = —0.039(7), for the point-split current 
and Z T Z S 1 = 0.185(7), for the local current ||. 
Also, the combination of ZtZ^ 1 shows no de- 
pendence on k. In fact, the latter is an 0(a) 
effect. An estimate of ZtZ^ 1 at the 1-loop or- 
der, using the point split current and (3 = 2.3 is 
Z T Zg X ee Z T \i-ioo P = -0.074 Q. Preliminary 



studies for ZtZ§ , using the local current defini- 
tion, are reported in |5^ , ^6| . Turning to Fig. |J, 
the expectation is that muZg 1 vanishes linearly 
when k — * k c . Also, a determination of k c by per- 
forming an extrapolation to zero gluino mass is 
obtained. The results is k c = 0.19750(38) for the 
point split current and k c — 0.19647(27) for the 
local one. These values can be compared with 
the previous determination from the first order 
phase transition, k c = (0.1955 ± 0.0005) @, see 
also Fig. [|. A study of the continuum limit of 
the lattice SUSY Wis, using Wilson fermions, at 
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Figure 3. ZtZ§ 1 as a function of 1/k with the op- 
erator insertion x(x) (filled diamonds) and Tq(x) 
(filled triangles). 
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Figure 4. am s Z~^ v as a function of 1/k with 
the operator insertion x(x) (diamonds) and Tq(x) 
(triangles). 



the 1-loop order is j55|,[56| . It is interesting to see 
whether an analytic calculation of the Zs and Zt 
renormalization factors is in agreement with the 
numerical ones. Compared with numerical simu- 
lations |33|, a lattice perturbative calculation of 
the SUSY Wis is more complicated. In order to 
do perturbation theory, we have to fix the gauge, 
which implies that new terms appear in the SUSY 
Wis: the gauge fixing term, the Faddeev-Popov 
term and a term coming from the gauge varia- 
tion of the involved operators J56fl , This is also 
discussed in J57|. Taking into account all these 
contributions we have |56| , 

(OVMx)) - 2m (O X (x)) + < J^r| E =o) - 



shell regime |56[ , give also results of the same or- 
der for Zt 



,05S GF . ,05S F p 



8e(x) 



Se(x) 



Se(x) 
| B=0 ) = (OX s (x)) 



In (55j|3(| the operator insertion studied (at the 
1-loop order) is O := A^(y) X b (z), which is a non- 
gauge invariant operator (a similar choice of op- 
erators as in the numerical case, for example x> 
would lead to a 2-loop calculation). In this case, 
Xs(x) mixes with operators of equal or lower 
dimension (3^,Q, as for the on-shell case, plus 
5^- ZiAi, where the additional operators A4 cor- 
respond to mixing with non-gauge invariant oper- 
ators. In [p5| , an on-shell 1-loop perturbative cal- 
culation (using an exact expression for Xs), gives 
ZrZg 1 = Zx\i-ioop = —0.074, using a point-split 
definition for the supercurrent, while preliminary 
studies using the local supercurrent, in the off- 



2.4. Ginsparg- Wilson (GW) fermions 

The GW relation (59) is defined to be 

j 5 D + D 75 = aD lb D , 



(21) 



where 75-D is an hermitian lattice Dirac operator. 
An explicit solution to the GW algebra, free of 
doubling species is reported in fl60j . It exhibits 
highlight chiralities properties plfl and locality 
p| . Either DWF (already discussed) or opera- 
tors obtained from a renormalization group (RG) 
blocking of the continuum one, the fixed-point 
(FP) operator [ p3[ , satisfy the GW relation. 

Theoretical applications of GW algebra to lat- 
tice supersymmetry are reported in p4|,|65] 66 
for the Wess-Zumino model fl68[ . In |34j, it is 
suggested that a lattice version of a perturba- 
tively finite theory preserves supersymmetry to 
all orders in perturbation theory, in the sense 
that the SUSY breaking terms induced by the 
failure of the Leibniz rule become irrelevant in 
the continuum limit. Differences between imple- 
menting GW fermions and Wilson fermions are 
also analyzed. In |35| a conflict between lattice 
chiral symmetry and the Majorana condition for 
a Yukawa- type coupling is pointed out. If one 
adopts GW operators, a precise analysis of super- 
symmetry and its breaking require a consistent 
formulation of Majorana fermions. 

In Ref. |6£j] a new formulation of SYM on the 
lattice with an exact fermionic symmetry is pre- 
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sented. First, it is considering the model in a 
fundamental lattice which is called the one-cell 
model and it is derived the preSUSY transforma- 
tions. Then, it is extended to the entire lattice. 
The lattice action has a peculiar form: a 2a trans- 
lational invariance, not the usual a ones, and it is 
called Ichimatsu lattice (similar to a chessboard) . 
At this conference, first non-perturbative results 
on an Ichimatsu lattice gauge theory is presented 
|70| , while in [[n] the study of the phase transi- 
tion is reported. Although this is an interesting 
formulation, one has to note that, because of the 
staggered fermion action, there is no exact bal- 
ance between fermionic and bosonic degrees of 
freedom. 



3. EXACT SUPERSYMMETRY 
THE LATTICE 



ON 



Improving lattice supersymmetry seems to be 
a difficult task for gauge theories. In fact, most 
SUSY theories, as for example N — 2 or N — 4 
SYM, contain scalar bosons which typically pro- 
duce SUSY violating relevant operators, which 
has to be fined tune away. For N = 1, as we have 
seen, only a fine tuning is needed in order to elimi- 
nate the mass term. Unfortunately, because there 
is no discrete version of supersymmetry which 
can be implemented to forbid scalar masses and 
unwanted relevant operators, it is desiderable to 
construct lattice structures which directly display 
at least a subset of exact supersymmetry in order 
to decrease the number of fine tuning to do. In 
Ref. 



72 



a similar approach is used, for the 2d 
Wess-Zumino model with extended N = 2 super- 
symmetry on the lattice. The model is then dis- 
cretized in a way that preserves exactly a subset 
of the continuum SUSY transformations. This is 
enough to guarantee that the full symmetry is re- 
stored without fine tuning in the continuum limit. 
Also numerical results using the HMC algorithm 
are presented Jt^] , and show the equality between 
fermions and boson masses and also the verifica- 
tion of the Wis to high precition. 

Another nice example of exact lattice super- 
symmetry is where a perfect SUSY action 
for the 2d and Ad Wess-Zumino model, free case, 
is presented. The perfect action is achieved in 



terms of block variable RG transformation, which 
maps a system from a fine lattice to a coarser lat- 
tice in a specific way, that keeps the partition 
function and all expectation values invariant. It 
is preserved invariance under a continuous SUSY 
type of field transformations in a local perfect lat- 
tice action, which contains also a remnant chiral 
symmetry. This perfect formulation also cures 
the problem related with the Leibniz rule on the 
lattice [22]] . In the perfect lattice formulation it 
is obtained the consistently blocked continuum 
translation operator. Therefore the algebra with 
the field variations closes. Ref. |6^| uses a similar 
approach. 

At this conference a talk by Kaplan |74|] , based 
on a recent remarkable paper |75| has been pre- 
sented. It is a new method for implementing 
supersymmetry on a spatial lattice for a variety 
of SYM theories in 3 + 1 dimensions, including 
N = 4, in a way that eliminates or reduces the 
problem of fine tuning. The formalism is pre- 
sented in the Minskowski space but a general- 
ization to Euclidean space is under way. The 
motivation is based on a recent work on decon- 
struction of SUSY theories (7(|. The spatial lat- 
tice is created by "orbifolding" . Starting from 
a "mother theory" , being a SUSY quantum me- 
chanics (QM) system with extended supersym- 
metry, a gauge group U(kN d ), and a global R 
symmetry group Gr, the spatial lattice is con- 
structed by orbifolding out a Zf^ factor from 
Gr x U(kN d ), which will result in a U(k) gauge 
symmetry, living on a <i-dimensional spatial lat- 
tice. This is called the "daughter theory" . Taking 
then the continuum limit of this "daughter the- 
ory", in some point in the classical moduli space 
of vacua, produces a higher dimensional quan- 
tum field theory with the original extended su- 
persymmetry restored (previoulsy broken by the 
orbifolding procedure), together with Poincare in- 
variance. The important point is that the lattice 
generated by orbifolding can retain some of the 
exact supersymmetries which facilitates the re- 
covery of the remaining ones in the continuum 
limit. This method also protects the renormaliz- 
ability of the theory. In a SUSY QM model 
with extended supersymmetry in + 1 dimensions 
is considered, and then the types of lattices that 
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can be obtained via orbifolding are shown. 

3.1. Other related topics 

At this conference, a nice example of a Hamil- 
tonian lattice version for the 2d Wess-Zumino 
model is presented by Campostrini [j77|. Previ- 
ous results in |7£| . Developing numerical simula- 
tions techniques using the Hamiltonian approach 
fn| turns out to be very advantageous (but con- 
tains also some disadvantages). Powerful many- 
body techniques are available: the Green Func- 
tion Monte Carlo (GFMC) algorithm |o). Since 
H is conserved, it is possible to preserve exactly 
a Id SUSY subalgebra of the continuum N = 1 
SUSY algebra. This is an advantage in compari- 
son with the standard lattice formulation and is it 
enough to guarantee the most important proper- 
ties of SUSY-like paring of positive energy states 
fr7| , frq| . Moreover, fermions are implemented di- 
rectly, and need not be integrated out. A disad- 
vantage is that fermions may lead to sign prob- 
lems, but at least in 1 + 1 dimensions, it can be 
by-passed. In |jT7| , simulations using GFMC are 
performed. The algorithm is very efficient in com- 
puting the ground state energy Eq (because it can 
distinguish between and 10~ 5 ), and therefore 
can be used to study the pattern of supersymme- 
try breaking, which corresponds to E > 0, while 
unbroken supersymmetry corresponds to E Q = 0. 
Focusing on V{(j)) = \(f> 2 + Ao, predictions for the 
strong coupling and the weak coupling regime are 
quite different and it is interesting to study both 
numerically and analytically and the crossover 
from the strong to the weak coupling Q. 

At this conference a related and very inter- 
esting method, namely, SUSY Discretized Light- 
Cone Quantization (SDLCQ) has been reported 
by Trittmann |pl[ . SDLCQ is a discrete, Hamil- 
tonian and manifestly SUSY approach to solving 
a quantum field theory. Some results (including 
a review) for the N = 1 SUSY Chern Simons the- 
ory, in 2d and 3d, including correlators and bound 
states are in |$2) . The theory is discretized by im- 
posing periodic boundary conditions on the boson 
and fermion field in terms of discrete momentum 
modes, k = nP/K, where K is a positive inte- 
ger that determines the resolution (typically of 
order 10 in the numerical calculations). The con- 



tinuum limit is reached for K — > oo. The low 
energy spectrum mass is determined using the 
fit M 2 = + 6(1/ K) and approximate BPS 
states which have non-zero masses are found. In 
f8lf , a numerical test of the Maldacena conjec- 
ture within the 10 — 15% is found. A schematic 
comparison of SDLCQ and lattice gauge theory 
is also presented. 

Ref. Q casts a Hamiltonian study of SYM 
QM, which faces the problem of the introduction 
of a cut-off that violates supersymmetry and the 
restoration of supersymmetry when it is taken 
away. It has been applyed to the Wess-Zumino 
QM and for the SYM QM, in 2d and id. There 
are no sign problems. Also, the complete spec- 
trum, Witten index and identification of SUSY 
multiplets has been determined Q . 



4. CONCLUSIONS 

A big effort has been made in order to describe 
supersymmetry on the lattice. Traditional Wil- 
son fermions have been used in realistic computa- 
tions with nice results. Improved chiral fermions 
results are starting too. New exciting ideas |75| 
are now waiting for numerical applications. It has 
also been shown that most properties of N = 1 
SYM are analogous to non-SUSY YM and QCD 
and can be tested on the lattice in order to under- 
stand a possible transition from SYM to YM. For 
this reason, placing real-world QCD in a wider 
variety of theories may help us to better under- 
stand it. 
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